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A REVIEW OF TURBULENCE MEASUREMENTS IN 
COMPRESSIBLE FLOW 

By 

V. A. Sandborn 

A survey of turbulence measurements in compressible flows is presented. 

The majority of turbulence measurements at super- and hypersonic speeds have 
been made for the zero pressure gradient, turbulent boundary layer. It was found 
that the nondimensional turbulent stress terms for the zero pressure gradient 
flow appear to agree closely with equivalent incomoressible measurements in 
the outer part of the boundary layer. The stress terms were nondimens ionalized 
by the wall value of shear stress and plotted versus the distance from the wall, 
nondimensionalized by the boundary-layer thickness. Indirect evaluation of 
the total shear stress distribution from mean velocity measurements for both 
super- and .hypersonic flows (zero pressure gradient , two 7 dimensional flows) 
indicate a near universal distribution. These total shear stress curves also 
agree very closely with measured incompressible shear stress distributions. 

Recent laser anemometer measurements of the turbulent Reynolds shear stress 
(puv) , reported by Johnson and Rose for a Mach number 2.9 flow, are in 
reasonable agreement with the expected total shear stress curve over the 
outer 60% of the boundary layer. Near the surface a systematic deviation 
between the laser measurements and expected shear stress was found. 

The laser anemometer data are only slightly higher than the incompressible 
results for the individual longitudinal and vertical turbulent velocity 
components over the complete boundary layer. These results might be taken to 
suggest, with some reservation, that the term (puv) may not be the only important 
term in the turbulent shear stress for supersonic boundary layers. 
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Hot-wire anemometer measurements of the longitudinal component of the 
turbulent velocity in super- and hypersonic flow were found to agree with the 
incompressible data in the outer region of the layers. Measurements of the 
vertical component of the turbulent velocity and the Reynolds shear stress 
with hot wires appear questionable. Mass flux and temperature fluctuation 
measurements show considerable scatter from one boundary layer to another. 

Neither mass flux or temperature data could be correlated with the degree of 
similarity observed for the turbulent stress terms. 

The spectral energy distribution of the super- and hypersonic flows are 
also shown to correspond to equivalent incompressible results. The spectral 
curves may be used to specify the frequency range of interest for a given 
boundary-layer measurement. Intermi ttency of the outer edge of the super- and 
hypersonic boundary layers is markedly different from that of the incompressible 
layer. . •••■*' 

Reported measurements of the free-streara mass flux intensity were found 
to increase as the square of the Mach number for both the super- and hypersonic 
wind tunnels. Very few measurements for subsonic compressible flow have 
been reported. A limited number of data points for the turbulent shear 
stress in a subsonic, compressible, pipe flow were found. 
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INTRODUCTION 

Techniques of dealing with turbulent shear flows are limited due to a lack 
of a rational theory of turbulent flow. Engineering calculations of flow about 
bodies still rely heavily on empirical information about the turbulence. Thus, 
it is necessary to experimentally evaluate turbulent shear flows, in order to 
upgrade and check the analytical methods. 

Measurements of turbulent properties in compressible flows have proven 
quite difficult. Although the basic concepts of measuring turbulence with 
hot-wire anemometers in compressible flows were evolved by Kovasznay, ref. 1 
and 2, in the early fifties, very few measurements were reported. In the late 
fifties a limited number of surveys by Morkovin and Phinney, ref. 3, and by 
„ Kistl.er, ref. .4, for supersonic boundary layers were, reported. Only within 
the last few years have further measurements of compressible turbulence appeared 
All the measurements are limited in content, so that it is difficult to access 
the accuracy of any one set of data. 

The present paper was undertaken as a comprehensive survey of compressible 
turbulence. In 1962 Morkovin, ref. 5, reviewed the effects of compressibility 
on turbulent flow and proposed limited guidelines as to expected effects. 
Morkovin showed that the effects of compressibility were expected to be somewhat 
passive. Based on the observations of Morkovin, together with recent observa- 
tions of large Reynolds number similarity of incompressible turbulent boundary 
layers, it was possible to examine in surprising detail a great amount of 
compressible turbulence data. The present survey is mainly limited to the 
simple flat plate boundary layer. However, it should serve as a guide to the 
evaluation of measuring techniques for studies of more complex flows. 



iv 


A 


e 


E 
E f 

f 

F(k) 

k 

L 

x 

M 

P 

P f 

Pr 

r 

Re 


R 


R 


mT 


R 


R 


OT 


T 

T » 


u 

U 

U 

T 


SYMBOLS 

supersonic tunnel geometry constant, eq. (8) 
hot wire anemometer voltage fluctuation 
total energy 

fluctuating part of the total energy 
frequency 

wave number spectral energy function, eq. (15) 
wave number 

macroscale of the longitudinal turbulent velocity comnonent in the 
x-direction, eq . (16) 

Mach number 

static pressure 

fluctuating part of the static pressure 

Prandtl number . * 

radial coordinate . 

Reynolds number 

Reynolds number based on the boundary-layer momentum thickness 

correlation between the mass flow and temperature fluctuations 

correlation of the longitudinal turbulent velocity components in the 
x-direction 

correlation between the entropy and vorticity fluctuations. 

turbulent heat flux 

temperature 

fluctuating part of the temperature 
x-direction turbulent velocity component 
x-direction mean velocity component 
shear stress velocity (/t^/p) 
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V 


v y-direction turbulent velocity component 

V y-direction mean velocity component 

w z-direction turbulent velocity component 

x longitudinal coordinate in direction of the mean flow 

y vertical coordinate in direction perpendicular to the mean flow and 

to the surface (for boundary layers) 

y* nondimensional vertical coordinate (yU_/v) 

z horizontal coordinate in direction perpendicular to the mean flow and 

parallel to the surface (for boundary layers) 

y intermi ttency factor (percent of time turbulent) 

5 boundary- layer thickness 

X microscale of the longitudinal turbulent velocity component 

p fluid density 

p 1 fluctuating part of the fluid density 

t shear stress 

x turbulent shear stress . 

t 

x wall shear stress 

w 

p coefficient of viscosity 

v kinematic viscosity 

Subscripts 

1 initial conditions 

e edge of shear flow 

f at the specific frequency f 

o total temperature 

r recovery temperature 

t turbulent part 


w 


wall value 



General Requirements for Turbulent Measurements 
in Compressible Flows 

Major interests in turbulent measurements arise from the need to evaluate 
flows about bodies, and in their wakes. These flows can be classified as 
turbulent shear flows. Other areas of interest include regions of turbulent 
mixing, as associated with problems such as jet engines and combustion processes. 
Also, since most aerodynamic studies are carried out in wind tunnels, a good 
deal of effort has been expended to evaluate the free-stream turbulence level 
of supersonic wind tunnels. Examination of the equations governing compressible, 
turbulent, flow will give an indication of the important turbulent quantities 
that need to be determined. 

Detailed derivations of the compressible equations are given by Van 
Driest, ref. 6, and by Schubauer and Tchen, ref. 7. For the present survey 
of compressible turbulence measurements it was found to be adequate to employ 
only the boundary- layer approximations of the complete equations. Following 
the analysis of Schubauer and Tchen, the steady-state boundary- layer equations 
may be written as: 

MOMENTUM 

\ C^UV J^t/W = r r ' 

CONTINUITY 

§ y {^ V ~^) =0 


(1) 

( 2 ) 

(3) 
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ENERGY 

k(*m * %(* &)=jy(fE- v+ i/i - ')M*f (4> 

Schubauer and Tchen argue that the turbulent shear stress, is just equal 

to ouv, which is equivalent to the incompressible results. Van Driest employs 
a slightly different approach wherein U p v remains lumped with the correspond- 
ing inertia term. The turbulent shear stress, for Van Driest f s analysis 

is given as 

q. S - (yOu/y — ~/5 Cw - Vy=>U u. v (5) 

It should be noted that evaluation of the total or turbulent shear stress from 
mean f low measuremen ts , assuming simila r ity ,' p r oduces a value equivalent to — 
that. of equation (5). Schubauer and Tchen neglect all but c uv. 

From the equations it is found that the following turbulence quantities 
should be specified. 

MOMENTUM U V ,0'V . sO'iX ^~ZCv 

■f ; s ) ' ' 

ENERGY (q W‘ - ~yO VET - V — yo’vE' 

Not all of these terms have been measured, nor can their magnitudes be estimated 
from the existing data. 

As is well known, equations (1) through (4) are not a closed set, due to 
the addition of the turbulence terms. Thus, the analytical evaluation must 
depend on some independent knowledge of the turbulence quantities. In particular, 
the shear stress and the heat flux terms are of major importance for engineering 
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predictions. For most boundary-layer flows it is found that the y-direction 
momentum equation, eq. (2), is independent of the x-direction equation, eq. (1) . 
Thus, equation (2) can be integrated directly to give 

c V 2 = (p - p ) (6) 

e 

If v- is small, then the variation of static pressure, p, across the shear 
layer can be neglected. 

For the boundary-layer evaluations that are considered later the vertical 
mass flux, o T v, is lumped directly with pV. Direct measurements of c r v 
do not appear to have been successful, so very little is known about its 
magnitude. Indirectly, o’v , could be evaluated from the continuity equation, 
eq. (3), by measuring the derivatives of the mean mass fluxes (for steady 
flow, |^=0) 

(7) 

The major interest in compressible turbulent boundary layers will be the 
evaluation of the shear stress. A number of measurements of o uv have been 
reported within the last year. Also, a large number of indirect evaluations of 
the shear stress from mean flow measurements have been reported for a wide 
range of super- and hypersonic boundary layers. For the case of zero pressure 
gradient, turbulent boundary layers, a surprisingly good agreement between 

the mean flow evaluations of turbulent shear stress was found in the course of 


the present survey. 



The basic equations of motion and energy are insensitive to much of the 
turbulent fluctuations. As a result, it has been necessary to seek higher 
order turbulent energy equations, ref. 8, in order to explore the " turbulent 
structure." These higher order turbulent energy equations have been employed 
in recent years as model equations to determine the turbulent shear stress, uv. 
Unfortunately, the assumptions necessary to deal with the higher order turbulent 
terms introduce a new set of uncertainties. Due to the difficulties in measuring 
turbulence in compressible flows, it is doubtful that experimental evaluation of 
the higher order turbulent energy equations is feasible. 

Although, not necessarily part of the solution of the equations governing 
turbulent shear flows, a large number of turbulent momentum and energy terms 
can be considered. 


MOMENTUM 


Turbulent Mass 

Shear Flow 


Triple 

Correlation 


p u 2 * V (pu) 1 p T uv 

p v ^ U (pv) 1 p T uw 

o ^ U(pw) f p 1 vw 

p uv 


p uw 
p vw 


ENERGY 


p v E 1 


V p f E f 


where u, v and w are the three turbulent velocity components. As previously 
noted, some of these terms are important in the equations, while others appear 
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as derivatives that can be neglected. For a large number of turbulent flows 
of importance a quasi- symmetry , either plane or cylindrical, can be obtained 
by proper choice of coordinates, so that the cross-velocity correlations, 
uw and vw (or uv depending on the coordinates chosen), can be neglected. 

These cross products will be zero if no shear exists in the z (or w) direction. 

The mass flux terms are of particular interest in evaluating the free-stream 
conditions of the compressible flow wind tunnels. Turbulent density and temperature 
fluctuations are also of interest in free-stream evaluations. In supersonic hot- 
wire anemometry evaluations the density fluctuations are directly coupled with 
the velocity fluctuations. Thus, it is found necessary to relate the independent 
temperature fluctuation measurements to the density. The relation between 
density and temperature fluctuations requires some knowledge as to the pressure 
fluctuations, since they are coupled through the equation of state. Thus, for 
most measurements the pressure fluctuations are assumed to be neglectable, and 
the density is related directly to the temperature. Independent measurements 
of the pressure fluctuations within the boundary layer are desirable, but 
experimental methods have not been available to measure them other than 
at the surface. A great deal of information has been obtained on the 
basic structure aspects of the turbulent fluctuations, such as: 

(a) energy-frequency content 

(b) probability distributions 

(c) correlations and related quantities such 


as scales and convective velocities. 
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REPORTED TURBULENCE MEASUREMENT'S 

A general review of existing measurements of turbulence in compressible flows 
has been made. The review both demonstrates what turbulent quantities can be 
measured, and the areas inhere information is available. At the present time 
it appears that very little information exists for the area of subsonic com- 
pressible flow. In some areas, such as zero pressure gradient, supersonic, 
turbulent boundary layers, a great deal of information is available. The present 
section is divided into four parts, according to the particular flow: 

(a) Free-stream Turbulence 

(b) Fully Developed Pipe Flow 

(c) Boundary Layers 

(d) Free Shear Flow 

(a) Free-stream Turbulence.- Historically, a good deal of effort was * 
expended in the mid-nineteen fifties to evaluate the free-stream turbulence 
level of supersonic wind tunnels. An early NACA attempt to compare the turbulence 
level of different supersonic wind tunnels consisted of measuring the transition 

location on a "standard" 10° cone, ref. 9. The transition was determined by 
measuring the recovery temperature distribution along the cone. Similar type 

studies have more recently, ref. 10, 11 and 12, been evaluated in the hypersonic 
facilities. Fischer and Wagner, ref. 10, obtained correlation curves of transition 
Reynolds number versus the free-stream disturbance levels for >_ 5 at the 
edge of the boundary layer. The transition Reynolds number was found to vary 
approximately as the inverse of the free-stream, mass flow fluctuation (or 
sound disturbance level) for a constant Mach number. An examination of the 
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limited measurements reported by Sandborn and Wisniewski, ref. 13, for the 
NACA, "standard" 10° cone in a small, 6- by 6-inch, Mach number 3 wind tunnel 
indicate a similar linear inverse relation. Note that it is necessary to 
separate the Mach number and turbulence level effects before different tunnels 
are compared. 

Figure 1 is a plot of the mass flow fluctuation intensity versus Mach 
number reported by a number of experimenters. For Mach numbers below 2, Laufer, 
ref. 14, and Morkovin, ref. 13, found that the velocity fluctuations upstream 
of the sonic throat affects the free-stream turbulence in the supersonic flow. 

At a Mach number of approximately 2.5 the upstream velocity fluctuations could 
no longer be related to the test section turbulence. For a Mach number 
range from 1.6 to 5, Laufer, ref. 16, found from hot-wire measurements that the 
correlation coefficient between mass flow and total temperature fluctuations in 
the free-stream had values of approximately -1. Secondly, it was known that, 
neither static temperature ? or velocity fluctuations by themselves were 
sufficient to produce the high turbulence levels indicated by the hot wire. 

Based on the analysis of Kovasznay, ref. 2, it was apparent that the only simple 
fluctuating field consistent with the measurements was a pure sound field. For 
a pure sound field the isentropic relations between pressure, density and 
temperature (and their fluctuations) are valid. Evidence for sound domination 
of the supersonic free-stream has been reported for a number of facilities. 

The sound fluctuations are related directly to the mass flow fluctuations 
measured with the hot-wire anemometer. The only deviations appear in the 
hypersonic flow facilities, where inlet flows must be heated to high temperatures 
For the hypersonic facilities, ref. 11 and 17, large temperature, as well as 


mass flow fluctuations are encountered. 
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The measurements of Laufer, ref. 16, Sandborn and Wisniewski, ref. 13, and 

Donaldson and Wallace, ref. 18, all show that the mass flow fluctuation intensity 

decreases with increasing Reynolds number at a fixed Mach number. The decrease 

-25 

is estimated to vary as Re # . It was pointed out by Laufer, ref. 16, that 

the wall boundary layers, which generate the freestream sound fluctuations, 
become thinner as the Reynolds number increases. The wall boundary-layer effect 
was dramatically demonstrated by Laufer, ref. 16, by reducing the tunnel Reynolds 
number to where the wall layers are laminar. The point shown on figure 1 at 

4 

M = 4.5 (Laufer, Re/in = 2.6x10 ) shows almost an order-of-magnitude reduction in 

tunnel mass-flow fluctuations. Although the data of figure 1 fail to agree in 

all cases there is evidence to suggest that the larger size supersonic wind 

tunnels for a given Mach number will have the lower turbulence levels (M > 2) . 

The ratio of f ree-s tream area to boundary-layer perimeter should be an important 

parameter. ' 

The increase of the mass-flow fluctuation intensity with Mach number is a 

result of the sound field being proportional to the fourth power of the Mach 

number, ref. 16. A curve fit of the data of Laufer for M > 2 shows that the 

2 

actual mass-flow fluctuations increase as M . As an approximate curve fit, 
the following relation was obtained 

( 8 ) 

where the constant A accounts for the variation due to tunnel size (and/or 
origin of the sound source). For Laufer's data a value A = 0.0095 was estimated. 
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Curves for A = 0,0095 and Re/in = 90,000 and 330,000, corresponding to 
Laufer’s measurements, are plotted on figure L A curve for A = 0.017 and 
Re/in = 90,000 is plotted to demonstrate the effect of A. This latter curve is 
an approximate fit of the Sandborn-Wisniewski data for the 6- by 6-inch tunnel. 

The measurements of Kistler, ref. 4 shown on figure 1, may not represent a true 
value for the wind tunnel, as they were taken at the edge of the wall boundary 
layer. Kistler reported that the free-stream levels were below the noise 
level of his measuring system. 

The very high Mach number results, shown as the insert on figure 1, do not 
vary according to equation (8). The particular facility is such that the 
Mach number increases as the Reynolds number is increased (see ref. 10). Thus, the 
higher Mach numbers also represent higher Reynolds numbers. Both the variation 
with Mach number and with Reynolds number are greater than would be predicted 
by equation (8). Note also that the mass flow intensities are much lower than 
the extrapolation of the supersonic curves, shown on figure 1, would predict. 

(b) Quasi-fully Developed Pipe Flow .- Fully developed pipe flow has 
proven of great value in the evaluation of turbulence measurement techniques 
for the incompressible case. It is possible that similar results can be obtained 
for compressible flows. For a compressible subsonic flow there is a progressive 
decrease in density along the pipe, so that the local velocity must increase 
with downstream distance. It is, however, found that a quasi- full-developed 
pipe flow may be obtained. Kjellstrom and Hedberg, ref. 19, consider the quasi- 
fully-developed flow momentum equation as 


-jrdlJ — _ dP _ 
d* dx 


±d_ 
r dr 


(tr) 


(9) 
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The stear stress variation with radius becomes 


f f _ r dP 

2 - c(A 



( 10 ) 


For incompressible flow the second integral in equation (10) is zero, so that 
the shear stress varies linearly across the pipe radius. Durst, Launder and 
Kjellstrom, ref. 20, have evaluated equation (10) for a typical empirical velocity 
distribution, They find a progressively increasing nonlinearlity due to the 
velocity term as the Mach number increases from 0.5 to 0.7. 

Figure 2 shows the results of an experimental study reported by Gibbings and 
Mikulla, ref 21, for subsonic, compressible, quasi-fully-develoned, D'ipe flow. 

Both the shear stress distribution experimentally evaluated from a form of 
equation (10) , and direct hot-wire measurements of puv were obtained. The 
nonlinear, "compressible" effect on the computed total shear stress was evident 
for a Mach number of 0.51, although it was small. The hot-wire measurements 
were found to agree with the evaluation of equation (10) over the center region 
of the pipe. Near the wall the value of puv drops to zero somewhat faster 
than might be expected. These measurements are the only subsonic, compressible, 
turbulence data found. 

(c) Boundary Layers.- The major interest of compressible turbulence mea- 
surement studies has been the supersonic boundary layer. A review of the litera- 
ture showed a number of studies of zero pressure gradient, supersonic, turbulent 
boundary layers has been made. These studies range in Mach number from 1.7 to 
9.4. Temperature and mass flow fluctuations have also been reported for M = 20. 
No measurements that contain compressible effects at subsonic Mach numbers were 
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found. Serafine, ref. 22, reports longitudinal turbulent intensity measurements 
for a M = 0.6 flow, however, all effects of compressibility were neglected. 

In the section on general requirements it was noted that the basic turbulent 
terms required for evaluation of the continuity, momentum and energy equations 
were: uv, p f v, p T u,p 1 uv,p 1 E f , vE ? and p T vE ! . These quantities have nroven very 

difficult to measure directly with hot-wire anemometers. Figure 3 is a summary 
plot of measurements of puv that have been reported for supersonic, zero 
pressure gradient, turbulent boundary layers. All of the data on figure 3 have 
been reported within the current year, so it would appear that major electronic 
difficulties are not expected. The shaded curve, shown on figure 3, represents 
a "best estimate" of the total shear stress (viscous plus turbulent, as given by 
eq. (5)) which will be discussed in more detail in the following paragraph. The 
measured values of puv were obtained with hot wire and laser anemometers. The 
measurements of Rose and Johnson, ref. 24, with a laser system do not require 
any assumption regarding the pressure fluctuations. If puv were the only term 
of importance in the compressible turbulent shear stress relation, as is the case 
for incompressible flow, the laser measurements should agree with the shaded curve 
shown on figure 3. The laser data agree approximately with the best estimate 
shear stress distribution over the outer 60 % of the layer. Near the surface 
the systematic deviation of the laser data from the expected results poses a 
question on the accuracy of the data. While it is difficult to access the 
accuracy, it will be shown that the evaluation of u 2 " and V 2 * from the same set 

rr 

of data are in reasonable agreement with the incompressible results. Both 
and v 7 were found to be slightly higher than the corresponding incompressible 
measurements over most of the boundary layer. 
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A second possibility for the deviation of the laser measurments (of cuv 
from the expected total shear stress curve) is that the other terms in the 
compressible shear stress, eq. (5), become important near the surface. The 
viscous term was estimated from the mean velocity measurements, and found to 
be two orders of magnitude smaller than the measured cuv term. The turbulent 
term Vp f u would be expected to be small near the surface, since V is small. 

A rough estimate based on a linear variation of V, with (V/ T J^ - 0.1, indicated 
this term is less than 5% of puv. Also Vp ' u has an opposite sign compared to 
cuv, so. it would act to increase rather than decrease the deviation near the 
surface. It appears that the deviations occur in the region where the turbulent 
fluctuation /v^ is of the same order of magnitude as the mean vertical velocity 
component, V. The term o’uv is generally suggested to be quite small, although 
there is no direct experimental information available on its magnitude. 

Although direct measures of the turbulent shear stress for supersonic 
boundary layers are not satisfactory, indirect evaluation from the equations of 
motion appear to give reasonably consistent results. For a zero pressure 
gradient, steady flow, equation (1) can be solved for the shear stress 

X ^(A |y -yBCTv- V/tt -yO^Tv ) - t* -hf *) ~ (s 5 Vvtydy (11) 

It is possible to make sufficient mean flow measurements to evaluate the mean flow 
derivatives in the integral term. However, for the evaluations that have appeared 
for equation (11) , an assumption of "similarity" was employed. The similarity 
introduced by Meier and Rotta, ref 25, assumes that the velocity and temperature, 
or density, profiles in the x-direction are self-similar. Exact similarity for 
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zero pressure gradient flow requires that the wall shear stress is constant 
with x, and the boundary-layer thickness vary linear with x. It is experi- 
mentally found that the incompressible, zero pressure gradient, boundary 
layer approaches the similarity requirements at "large" Reynolds numbers, 

Zoric, ref. 26. 

Figure 4 is a summary plot of the total shear stress distribution) com- 
puted from similarity forms of equation (11) . The data shown in the main 
part of figure 4 are for measurements along flat plates, conical models and 
on two-dimensional tunnel walls. A wide range of both Mach number and 
Reynolds number are represented in the data used to compute the distributions 
shown on figure 4. There are surpisingly only small variations in the diverse 
number of measurements. The data shown are mainly for near adiabatic flow 
conditions, although heat transfer is also represented. The data of Horstman 
and Owen, ref. 27, Danberg, ref.' 2*8, (evaluated by Bushnell and ’Morris, 
ref 29) and Samuels, Peterson and Adcock, ref 30, (also evaluated by Bushnell 
and Morris) are all for wall to freestream temperature ratios of the order 
of 0.45 to 0.48. Of the data found in the literature only a set of three 
profiles reported by Rochelle, ref. 31, for M = 2.029, 2.480 and M = 4.975 
are not included on figure 4. The profile evaluated by Rochelle for M = 2.480 
falls approximately within the limits of those shown on figure 4 (it appears to 
indicate a decreasing pressure gradient), but the other profiles fall well 
below the present results. Apparently some question arises in the technique of 
calculation used in this early evaluation. 

A "best estimate" for the shear distribution for two-dimensional zero 
pressure gradient, supersonic, turbulent boundary layers was constructed 
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from the curves of figure 4. The "best estimate" distribution is shown on 
figures 3, 4 (insert) and 5 as the shaded regions. From the equation of 
motion, eq. (1) the boundary condition for the shear stress at the wall is 


ov 

" y=o 


£2.1 
dx ' 


y=0 


( 12 ) 


For a zero pressure gradient the shear distribution must approach the 


surface with zero slope. Of the distributions shown on figure 4, the pro- 

3 

files of Horstman and Owen, ref. 27 (for R = 9.7x10 only) and the one 

d 

from Danberg, ref 28, (evaluated by Bushnell and Morris, ref. 29) appear 

to indicate a negative slope for -~*J , Thus, in drawing the best estimate 

y=0 


distribution the region near the wall does not include these two dis tributions . 

The scatter between the profiles in the outer region of the boundary 
layer may represent an experimental uncertainty in defining the boundary- 
layer edge. The boundary-layer edge .is a statistical property, which is 
not well defined in the mean-flow measurements. Incompressible measurements, 
such as those of Klebanoff, ref. 32, for zero-pressure-gradient boundary 
layers, show that the turbulent shear stress goes to zero at values slightly 
greater than the statistical boundary-layer edge. At the boundarv-laver 
edge the incompressible turbulent shear stress was approximately two orders 
of magnitude less than the wall shear stress, as shown on figure 5. It is 
reasonable to assume that the supersonic shear stress should also approach 
closely to zero at y/6 = 1. 

Several shear stress distributions for boundary layers in nozzle- type flows 
are shown on the insert of figure 4. It is not at present possible to determine 
the validity of the similarity concept in the evaluation of these flows. 
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It is, however, doubtful that the large differences between axi symmetric and 

two-dimensional flows, shown on figure 4, are to be expected. 

Between the different sets of measurements represented on figure 4, it 

was not possible to identify consistent Mach number or Reynolds number 

variations. The systematic variation of Mach number at constant Reynolds 

number, reported by Meier and Rotta, ref. 25, appears to show a small 

systematic increase in the values of t/t^ with increasing Mach number for 

most of the boundary layer. With perhaps the exception of the M = 4.5 data, 

the variation with Mach number was quite small. Systematic variation of 

Reynolds number reported by Horstman and Owen, ref. 27, and also Meier 

and Rotta indicate a slight decrease in the values of t/t with increasing 

w 

Reynolds number for most of the boundary layer. The possibility exists 
that, since all the evaluations employ the similarity approximations, the 
resulting distribution must be locked to a similarity .shape. Obviously, 

the fixing of the distribution to T / T w = 1 at y/5 = 0 and ~ 0 at 

y/5 = 1, together with a zero slope at y/6 * 0, makes large variations 
unlikely. 

For subsonic flow the variation of the zero pressure gradient, turbulent 

shear distribution from the near-similarity condition to a flow with an order 

of magnitude smaller Reynolds number is within the accuracy of the measurements. 

Figure 5 compares the faired incompressible measurements of turbulent shear 

3 

stress reported by Zoric, ref. 26, for R^ = 42x10 with Klebanoff f s data, 

3 

ref. 32, for R = 7.75x10 . The flow evaluated by Zoric was very close to 

9 

the similarity type that is presently assumed in the supersonic shear stress 
evaluations. The comparison of the incompressible measurements appears to 
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support the use of the similarity assumptions to evaluate the turbulent 
shear stress. The incompressible measurements might well represent the lower 
limit of the "best estimate" distribution. 

Maise and McDonald, ref. 33, have derived shear stress distributions 
based on an empirical correlation of measured boundary layer, mean velocity 
profiles. Their predicted shear distribution for M = 5 is shown on figure 
5. The empirical distribution is higher than the experimental evaluations 
over most of the layer. A similar over-prediction of the incompressible 
shear distribution is also noted in their results. The analysis does, however, 
predicte a small increase in the shear stress as Mach number increases, and 
also a small decrease in shear stress as Reynolds number increases. 

In summary it may be proposed that measurements of the turbulent shear 
stress, (-Quv -o ? uv - Vp f u), for supersonic, zero pressure gradient, two- 
dimensional, adiabatic,, turbulent boundary layers should fall within the 
"best estimate" curve, shown on figure 3. The exception, of course, is the 
region very close to the surface where the turbulent shear stress must drop 
to zero. The zero pressure gradient flow should serve as a check of the 
turbulence measuring technique. 

Near the surface the viscous shear stress can be computed from the mean 
velocity measurements, so that the turbulent shear stress drop-off can be 
predicted. For incompressible flow of air the region of turbulent shear 
stress drop-off occurs extremely close to the surface, which makes actual 
measurements difficult. The supersonic turbulent boundary layers tend to 
have thicker viscous sublayers, which may help to make measurements more 
accessible. The sublayer becomes increasingly thicker as the Mach number 
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increases. In the past, hot-wire measurements of turbulent shear stress near 

the surface in incompressible flows produce values that are higher than 

expected, Sandborn, ref, 34, p. 294. The measurement difficulty was found 

to be due to the yawed hot wire being in a large turbulent velocity gradient, 

so that nonlinear averaging across the wire occurs. 

The present analysis has demonstrated that an ’’outer region” boundary- 

layer similarity can be used to correlate the total shear stress distribution 

over the complete layer. It is also generally accepted that an ’’inner region” 

similarity may be used near the surface. For subsonic flow the velocity near 

the surface is found to scale as the ”friction/or shear velocity,” U_(^ J • 

^ w 

Thus, for the region where the viscous shear stress is important, the turbulent 
shear stress, as well as the other turbulent stress terms, might be expected 
to be a universal function of the characteristic velocity, U T and length 

yU i • 

y*( = — . Requirements for the existence of similarity for the^ outer region 
(U = constant, 36 /3x = constant) are more restricted than for the inner 

T 

region similarity (3U^/3x = constant) (see ref. 26). For super- and hypersonic 
flows the ’’inner region” similarity can be expected to become increasingly 
important. The extent of the viscous region increases rapidly with increasing 
Mach number. For a Mach number of 5 the sublayer, defined as y* 11, is 
roughly ten times larger than for the same Reynolds number at Mach number zero. 

It is found that the mean velocity measurements are correlated by the 
inner similarity coordinates for all Mach numbers within the sublayer (y* < 11) . 
Proper inclusions of the density variation, such as that given by Van Driest, 
ref. 6, are found to correlate the mean velocity data out to values of y* 
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greater than 200. For subsonic flow the turbulent velocity components and shear 
stress are found to be correlated by the inner similarity coordinates only in 
the viscous sublayer. Figure 6 shows a set of turbulent shear stress distributions 
obtained by Laufer, ref. 35, from mean velocity and pressure measurements in a 
fully developed pipe flow. Note that the inner region similarity for this ideal 
flow (U_ » constant) only extends out to y* - 20. A set of hot-wire measurements 
of the turbulent shear stress reported by Tieleman, ref. 36, for the large 
Reynolds number, incompressible boundarv layer is also shown on figure 6. It is 
difficult to obtain direct measurements of the turbulent shear stress in the 
wall region. The hot-wire measurements of Tieleman required up to 20% corrections 
due to the nonlinear effects of turbulent velocity gradients along the wire. 

The measurements are in agreement with the results of Laufer. 

Shear stress measurements in supersonic flows do not reach into the viscous 
sublayer. The turbulent shear stress distributions determined from mean velocity 
measurements by Horstman and Owen, ref. 27, were employed to evaluate the 
compressible effects on the turbulent inner region similarity. These calculations 
were made assuming the mean flow is similar in the outer region. The insert 
on figure 6 shows the comparison of four sets of data taken along a cylindrical 
model in a Mach number 7.2 flow. For the wall coordinates used in the insert 
the similarity does not exist beyond y* - 13. The main figure 6 is a renlot 

T T 

of the data in terms of — = , where - p is akin to a compressible 

T w p (y)ir p w p w 

scaling parameter originally suggested by Morkovin, ref. 5. The Van Driest 
correction is somewhat more complex in that an integral correction for density 
is required. No correction is required by the y* coordinate. The density 
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correction correlates the data over the region 10 < y* < 40, which is somewhat 
better than found in Laufer’s pipe flow. However, due to the limited region 
of the inner similarity, it was not possible to compare existing compressible 
turbulence measurements in these coordinates. The disagreement in magnitude 
between the results of Horstman and Owen and incompressible results may be due 
to uncertainties, such as the extrapolation of the outer region similarity to 
the wall. 

The present review demonstrates that the general magnitude of the shear 
stress distribution across a supersonic, adiabatic, zero pressure gradient, 
turbulent boundary layer can be estimated to within approximately +15%. The 
general observation of Morkovin, ref. 5, that compressible effects on the 
turbulence are "passive" appear valid for the shear stress distribution. It is 
possible that the incompressible, zero pressure gradient, shear stress distri- 
bution is an adequate model for the -super- and hypersonic flows, at least up to 
Mach number 7. 

For other than the "ideal" flat plate case information on the turbulent 
shear stress is quite limited. No shear stress information for highly nonadiabatic 
cases were found. Maise and McDonald, ref 33, demonstrated the nonsimilar nature 
of the mean velocity profiles for several cases, which suggest major variations 
may occur. A limited amount of information is available for pressure gradient 
boundary layers. Rose and Johnson, ref. 24, and ref. 37, have reported hot-wire 
and laser anemometer measurements of puv for shock-wave- induced, positive pressure 
gradient flows. Sturek, ref. 38, gives calculations of the shear stress dis- 
tributions (from evaluation of the equation of motion) for an isentropic-ramp 
compression flow. In subsonic flow a positive pressure gradient, which is 
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referred to as an adverse pressure gradient, normally produces a decrease in 
wall shear stress and the turbulent boundary layer approaches separation. 

Although, special adverse pressure gradients have been produced with constant 
or slightly increasing wall shear stress. For supersonic flow a positive 
pressure gradient can produce either (or both) boundary-layer separation and/or 
a '’compression" of the layer with large increases in the shear stress both at 
the wall and in the layer. The "compression” type boundary layer decreases in 
thickness and developes a maximum in the total shear stress that can be many 
times greater than the wall value. It is difficult to conceive of a subsonic 
flow that is equivalent to the supersonic "compression" layer. Supersonic boundary 
layers experience the adverse pressure gradient effect of decreasing wall shear, 
which may result in flow separation, slightly upstream or in the region of a 
shock-wave-interaction or a compression turn. However, once through the shock 
w^ve or on an isent ropic-ramp the layer is subjected to the positive pressure 
gradient compression effect. 

Figure 7 is a plot of the turbulent shear term, puv, referenced to the wall 
shear, measured by Rose, ref. 37,* and by Rose and Johnson, ref. 24. The data 
of Rose is divided into three regions (figures 7a), b) , and c)), the approach 
region, the shock wave interaction region, and the compression region. The 
wall shear stress increases with x-distance through the interaction and compression 
regions. Although the approach data for puv do not agree with the expected 
values of the total stress, the general increase in value through and downstream 
of the interaction are likely to be characteristic of a "compression" layer. 

The laser measurements of Rose and Johnson, ref. 24, (figure 7d)) show the 


*[Rose plots the term (Puv 4- Up'v) as the turbulent shear stress in his report.] 
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characteristic large increase in the total shear stress in the compression 
region. The pressure distribution downstream of the shock-wave interaction, 
reported by Rose and Johnson, is nearly constant, so the boundary layer at 
x = 9,375 should just be starting to recover to the zero pressure gradient case. 

Figure 8 shows the total shear stress distributions evaluated from mean 
flow measurements by Sturek, ref. 38, for an isentropic-ramp -induced adverse 
pressure gradient flow. Although it was not clear from the paper it appears that 
the evaluations are not based on any similarity assumptions (some reservations were 
expressed on the assumption of a constant static pressure through the layer).* 

Both the analysis of Sturek and the measurements of Rose and Rose and Johnson 
indicate large increases in the magnitude of the turbulent shear stress at some 
distance away from the wall. The compression flow produces the increase in the 
wall shear stress, and also equation (12) requires that the shear stress increase 
away from the wall. In contrast, for incompressible, adverse pressure ■ gradient 
flow a decrease in magnitude of both the wall shear stress and the maximum 
turbulent shear stress is normally observed. The incompressible results of 
Sandborn and Slogar, ref. 39, are included as an insert on figure 8 (where both 
profiles are referenced to the initial wall shear stress value). Thus, the 
supersonic "compression" flow is quite different from an incompressible, adverse 
pressure gradient flow. The compression effect proves to be a good "generator" 
of turbulent shear stress. 

It is expected that a supersonic turbulent boundary layer approaching separation 
would have a shear stress distribution that varies much as shown in the insert 


*The large variation of t with distance downstream may also raise the question 
of the validity of the boundary-layer assumptions. 
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on figure 8. No turbulence data on a supersonic, separating, turbulent boundary 
layer was available. The actual variation from the zero pressure gradient 
shear stress distribution to the separation distribution obviously must occur 
over a very short distance. 

Evaluation of the turbulent heat flux term, = (-pvE 1 - Vc ! E ! -c 1 vE 1 ) (13) 

can be obtained in much che same manner as the shear stress term. For the 


insulated flat plate flow Meier and Rotta, ref. 25, have determined the turbulent 
heat flux from the energy equation, eq. (4). Similarity assumptions are also 
required for the heat flux evaluations. Figure 9(a) is a plot of the turbulent 
heat flux evaluated by Meier and Rotta as functions both of Mach and Reynolds 
number. The turbulent heat flux is found to increase with increasing Mach number , 
and decrease with increasing Reynolds number . (For the different Mach numbers 
shown, the Reynolds number (Re = *1.9x10 5 ) is approximately constant. 

Figure 9b) .shows the evaluation by Horstman . and Owen, ref. 27, for the boundary 
layer along a cylinder with heat transfer, T^/T q = 0.46. The Reynolds number 
effect is also seen in the heat-transfer case (i.e.» decreasing q^ with increas- 
ing Reynolds number). Positive values of q^_ (not shown on figure 9b)) are 
obtained for the heat transfer case near the surface, however, the region near 


the surface is uncertain due to the similarity assumptions, 
portion of the boundary layer the "turbulent" Prandtl number. 


Over the major 

T t (3E/ay) 

Pr ” q t (Sv/3y) 



was approximately 0.9 to 0.8 for both the adiabatic and the heat- trans fer case. 


Near the surface, where the similarity assumptions are least valid, the turbulent 


Prandtl number increases to values greater than one. Near the outer edge of 


the layers the turbulent Prandtl number drops to zero. Direct measurements of 


the turbulent heat flux terms have not been reported. 
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For the evaluation of hot-wire anemometer measurements in supersonic flow, 
it is usually found necessary to assume that the pressure fluctuations are zero. 

The assumption of zero pressure fluctuations together with the equation of state 
leads to the requirement that (see for example, Rose, ref. 37) 

:v'E f =-Ep ? v (14) 

For this case the turbulent heat flux is just equal to Vp'E 1 . Rose, ref. 37, 
reports measurements of the terms vT 1 and p T v, which are of experimental necessity 
consistent with equation (14) . 


The longitudinal turbulent velocity, has been measured by a number 

of experimenters. Figure 10 is a plot of the longitudinal turbulent stress, pu , 
referenced to ^ the local wall shear stress . (square root), versus y/q. *Since the 
total shear stress was found to be nearly similar, when referenced to the wall 
shear and the boundary-layer thickness (figure 4); it might also be expected 
that each of the terms in the turbulent stress tensor will show an equivalent 
similarity in the outer region. The incompressible measurements of Zoric, ref. 
26, demonstrated that all the terms in the turbulent stress tensor obtained 
similarity forms in the outer region. For the incompressible, high Reynolds 
number case evaluated by Zoric the similarity appears to hold for values of 
y/o > 0.05. The approach to similarity of the individual turbulent velocity 
components was found to be slower than that observed for the mean velocity or 
the total shear stress. The deviations of the high Reynolds number, similarity 
measurements of Zoric, ref. 26, and the low Reynolds number data of Klebanoff, 
ref. 32, in the outer region are within the experimental accuracy of the measure- 


ments. 


was originally proposed by Morkovin, 


/u 2 o (y) _ pu z 

The coordinate U . p lx 

T v W W 

ref. 5, in 1962 totally independent of the present similarity concept. Morkovin 

viewed j 0 ^ as the compressible scaling parameter. Thus, figure 10a) is 

v p 

w 

equivalent to figure 3 of Morkovin T s paper, except recent information has been 
included. 

All of the supersonic measurements shown on figure 10 appear to agree closely 
with the incompressible measurements in the outer part of the layer. Only the 
data of Johnson and Rose, ref. 23 and 24, appear to agree in the region close 
to the surface. The fact that the hot wire and laser measurements of Johnson 
and Rose agree with each other was taken as justification that pressure fluctua- 
tions could be neglected in evaluating the hot-wire output. The measurements 
of Kistler, ref. 4, and Rose, ref. 37, in the inner part of the layer show 
more of an inconsistency with the .data of Johnson and Rose than would be 
attributed to either Mach or Reynolds • number variations. 

The measurements of Owen and Horstman, ref. 40, and Laderman and Demetriades, 

ref. 17, shown on figure 10b) are for cooled wall cases of T /T =0.46 and 

e 

0.38 respectively. The stabilizing effect of the cooling can act to reduce the 
turbulence. Thus, it is possible that the heat transfer produces the large 
deviations from similarity seen in figure 10b) *. Pressure fluctuations may also 


*[Note that the total shear stress distributions evaluated by Horstman and Owen, 

ref. 27, and shown on figure 4, are for approximately the same flow conditions 

3 

as their data of figure 10b). The shear stress profile (R^ * 9.7x10 ) of figure 
4 is at the same station, x = 225 cm, as the turbulence measurements (circles) of 
figure 10b). The total shear stress profiles were computed assuming similarity 
exists. As noted, this profile does not give the expected zero slope at the surface, 
so that the assumption of similarity near the surface for the hypersonic flow 
of Horstman and Owen appear questionable.] 
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be a problem in the hypersonic flows, although Laderman and Demetraides have 
attempted to correct for the pressure fluctuations. 

Measurements of Rose, ref. 37. and Rose and Johnson, ref. 24, of the longi- 
tudinal turbulent stress variation through a shock-wave-boundary-layer interaction 
are shown in figure 11. The measurements in the approach region of the flows 
are included on figure 10a). The approach region longitudinal turbulent stress 
was in good agreement with the expected similarity. The effect of the shock-wave 
interaction, as seen from figure 11, is to increase the longitudinal turbulence 
in the inner two-thirds of the layer. The absolute magnitude of the longitudinal 
turbulent shear actually increases with x throughout the flow, since the value 

of t increases with x. The distributions maintain a definite similarity in 

w 

the outer region of the layer Any deviation of the measurements in the outer 
region* from the similarity profile was within the uncertainty in defining the 
boundary-layer thickness. The boundary- layer thickness had to be determined from, 
temperature measurements, since the shock wave masks the velocity variation. 

Figure 12 is a summary plot of reported measurements of the vertical turbulent 
stress, /ov 7 , nondimensionalized by the wall shear stress. Problems inherent 
in the measurement of the turbulent shear stress are also present in the measure- 
ment of & . The incompressible measurements also show large* disagreement than 
previously noted for the total shear stress or the longitudinal turbulent stress. 
The large Reynolds number boundary-layer measurements of Zoric, ref. 26, do 
not show the slow drop-off in /v^, as the wall was approached. Detailed measure- 
ments of the region very close to the surface in these large boundary layers were 
reported by Tieleman, ref. 36. The data of Tieleman, also included on figure 12, 
(note that the measurements were corrected for turbulence gradient effects on 
yawed wires) indicate a maximum value of /v^ occurs at a value of y/6 - 0.025, 
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which corresponds to a physical distance of 2,06 cm from the surface. The maximum 
value of / v^~ for Klebanoff f s data occurs at a value of y/5 - 0.18, or y - 1.37 cm. 

The laser anemometer measurements of Johnson and Rose, ref. 23, for supersonic 
flow were found to agree reasonably well with the large Reynolds, incompressible 
measurements of Zoric and Tieleman, as shown on figure 12. The measurements of 
Rose and Johnson are perhaps slightly higher than might be expected. The agree- 
ment for both the u and v components was taken as a strong indication that 
the laser measurements were reasonable, although it is not established that the 
supersonic results should definitely agree with the incompressible data. As noted 
previously, these results would seem to suggest that the laser measurements of 
ouv are also valid. 

Laser measurements of /:r downstream of a shock-wave-boundary-layer 
interaction were reported by Rose and Johnson,- ref. 24, and are shown in figure 
13. Their data- show no appreciable effect due to the interaction. Hot-wire - 
anemometry measurements through shock-wave-boundary-layer interactions are reported 
by Rose, ref. 37, and Rose and Johnson, ref. 24, However, it appears that questions 
still exist in the evaluation of hot-wire data for the measure of "v 2 ", as can be 
seen from figure 12. 

The basic measurement with the hot-wire anemometer in compressible flow is 
the mass flow fluctuation, (pu) 1 , since the heat transfer from the wire is 
equally sensitive to the density and the velocity. Thus, one would assume that 
the mass flow fluctuations are more accurately measured than any of the turbulent 
stress tensor terms. Figure 14 is a summary plot of the faired profiles of the 
mass flow fluctuations, referenced to the local mean mass flow, versus y/6 for 
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a number of supersonic and typersonic turbulent boundary layers. The data of 
Kistler, ref. 4, for one particular wind tunnel shows a systematic increase in 
the local mass flow fluctuations with increasing Mach number. The variation of 
the mass flow fluctuations from one boundry layer to another is too great to 
suggest any systematic effects. 

The mass flow fluctuations were reported by Rose, ref. 37, for the shock- 
wave-boundary-layer interaction. As might be expected, there was an increase 
in (pu) * through the compression region. Possibly the high levels of 
(cu)'/pU for the approach region, shown on figure 14, may be due to the nonuniform 
compression effects in the upstream nozzle flow. Values of c'u and o'w are 
also given by Rose. The values of p'w were nearly zero throughout the complete 
flow region. The values of p'v are roughly one half of p f u in the approach 
region and increase to about the same values as p ? u in the interaction and 
compression regions. ' • 

Temperature fluctuations have also been reported for most of the boundary 
layers where turbulence measurements were made. The resistance thermometer, 
which is an unheated hot wire, was employed for temperature measurements. In 
principle total temperature f .luc tuations can be measured directly without 
interaction with density, pressure, or velocity fluctuations. Unfortunately, 
difficulties are encountered, in that the thermometer supports may not be at the 
same mean temperature as the sensor. The heat- transfer error due to the supports 
can be large, particularly at the high Mach numbers and/or low density flows. 

Figures 15a) and b) are summary plots of the reported total and static 
temperature fluctuations. It is not readily apparent how to correlate the 
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different measurements. Morkovin, ref. 5, suggested the total temperature 
fluctuations be nondimensionalized by dividing by the difference between free- 
stream static and total temperature. The present correlation employs the 
recovery temperature rather than the total temperature, as a possible improvement 
in the inclusion of hypersonic data. The hypersonic data available, Owen and 
Horstman, ref. 40, and Laderman and Demetriades, ref. 17, include effects of 
heat transfer. It is not obvious that temperature fluctuations in the presence 
of heat transfer would be expected to agree with those in adiabatic flow. It 
is surprising that the heat- transfer flows do not indicate larger temperature 
fluctuations than the adiabatic flows. 

The static temperature fluctuation intensities are weighted bv a repre- 
sentative mean static temperature difference across the boundary layer. This 
weighting is similar to that employed by Kister, ref. 4, except the recovery 

temperature . is used (figure 15b) rather than the total temperature. These 

. * 

coordinates are based on the concept that the temperature fluctuation level should 
be proportional to the mean temperature gradient across the boundary layer, which 
was proposed originally by Kovasznay, ref. 2. With the exception of Rose's 
measurements, the supersonic boundary layers appear to be reasonably well correlated 
with one another. The hypersonic, heat transfer, boundary layer data agree 
roughly with the adiabatic static temperature variation. 

Figure 16 is a summary plot of measurements of the mass flux- temperature 
correlation. The correlations appear to vary greatly from one boundary layer 
to another. The best agreement is for the cooled wall, hypersonic, boundary 
layers, which appear to suggest a near constant value o f R mT 


- 0.6 over most 
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of the layer. The supersonic data reported by Rose, ref. 41, give correlations 
of approximately one. Evaluation of the correlation uT 1 , neglecting pressure 
fluctuations, was also reported by most of the experimenters. 

The major missing measurements in compressible turbulent boundary layers are 

the pressure fluctuations. As noted, it was necessary to make some direct 

assumption regarding the magnitude of the pressure fluctuations, in-order to 

evaluate the hot-wire measurements . Most of the evaluation techniques simply 

neglect the pressure fluctuations. Laderman and Demetriades, ref. 17, employed 

an assumption that the entropy and vorticitv were out of phase, but perfectly 

correlated (R^ = -1). Also they assumed that the sound mode was not correlated 

with either the entropy or vorticity modes. These assumptions, together with 

their data give pressure fluctuation intensities that are slightly larger than 

the corresponding density fluctuation intensities. In figure 10b) the measure- 
“T* 

ments of / - — evaluated by Laderman and Demetriades, ref. 17, with the 

* T * ' 

• W ‘ 

assumed corrections for pressure fluctuations are compared with measurements of 
Owen and Horstman, ref. 40, where the pressure fluctuations were neglected. The 
corrections do not appear to show a measurable effect at least in the middle 
region of the layer. The agreement of the hypersonic data with the incompressible, 
similarity profile of Zoric would appear to justify the neglect of pressure 
fluctuations in the outer region. 

While direct measure of pressure fluctuations in the boundary layer, other 
than at the surface, have not proven feasible, it is possible that independent 
measurements of density fluctuations can be made. Recent developments in optical 
schlieren techniques, Wilson and Damkevala, ref. 42, and interferometry techniques, 
Wehrmann, ref. 43, have shown the feasibility of measuring density fluctuations. 
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For low density flows, Wallace, ref. 44, has used an electron beam to obtain 
density fluctuations. Harvey and Bushnell, ref. 45, show the variation of 
the density fluctuations obtained from Wallace’ s measurements. These measure- 
ments were taken in a hypersonic shock tunnel, with a very large temperature 
difference across the flow. The measured density fluctuations are much lower 
than would be indicated by the temperature evaluations shown on figure 15b) . 

Measurements of the detailed structure of the turbulent fluctuations are 
also available. In order to evaluate the frequency requirements of the 
instrumentation used for super- and hypersonic turbulence measurements, the 
energy spectrum is necessary. Some form of the mass flux or longitudinal 
turbulent velocity spectra is reported for most studies, unfortunately, it 
is quite difficult to compare the many different forms of the spectra that have 
been presented. Theoretical considerations suggest a spectra form referenced 
to a characteristic wave number and dissipative scale be used. For large 
Reynolds numbers, it is found that both atmospheric and wind tunnel turbulence 
spectra reduce to a universal curve for the higher wave numbers, when these 
scaling parameters are employed, Sandborn and Marshall, ref. 46. At the low 
wave number or low frequency end of the spectra the physical scale of the 
local flow will determine the individual spectra. The required dissipation 
scales are not available for the high-speed flows studied, so the evaluation 
of the universal spectra is not attempted. As a rough approximation, it 
might be assumed that the flow scales for the high speed tests are about the 
same. Thus, the comparison of the measured spectra in wave number coordinates 
might be a logical choice. Detailed measurements of the spectral variation 
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through a large turbulent boundary layer for incompressible flow was reported 
by Tieleman, ref. 36. For the incompressible flow a family of systematic 
developing spectra were found near the surface, (Insert, figure 17). Once 
the outer region, where similarity of the mean flow exists, was reached the wave 
number spectra are all identical and do not vary with the y-distance. 

Figure 17 is a plot of the wave number spectrum; 

where 

03 

/ F(k) dk i 1 (15) 

and 

2TTf 

k = wave number = — — ; where f is the frequency 

in hertz, and U is the local 
mean velocitv in ft/sec 

u u f 

F(k) = x — — ; where u^ is the mean square 

u turbulent energy per unit hertz 

at the frequency, f, and 

u • is. .the total mean. square 
turbulent energy for all # frequencies 

The data of Horstman and Owen (not previously reported, but for the same flow 
as reported in reference 40) are for the mass flux, (pu) ! . The other spectra are 
for the longitudinal turbulent velocity, u 2 ". With perhaps the exception of 
the outer edge spectrum of Horstman and Owen (y = 3.36 cm and y/5 = .96), the 
spectra appear to vary directly as the distance from the surface. 

Figure 17 should serve as a means of determining the expected frequency 
range of interest in boundary-layer measurements. Once the mean velocity, 
as a function of y-distance, is known, the frequency corresponding to any 
given wave number can be computed. The lower bound on wave number might be 
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represented by the curve of Tieleman for y = 17.8 cm. The upper limit on 
wave number most certainly will fall within the limits of the curve given bv 
Tieleman for y = .0025 cm. This curve at y = .0025 cm is approximately the 
same as the data of Morkovin, shown on figure 17, for the lower wave numbers. 

Note the most severe frequency requirements occur when very high velocities 
(thin boundary layers) are measured very close to the wall. 

Spectra for temperature fluctuations are also available. The frequency 
content of the temperature fluctuations appear similar to the velocity fluctua- 
tions. Laderman and Demetriades, ref. 17, find that their temperature and 

velocitv spectra are nearly the same when plotted as up/ u z , - N and 

f (f=o) 

fL 

Ti 7 - /T' j ^ x versus - 77 - . The scale length L corresponds to Taylor's 
f (r=o) U 0 x 

macroscale of turbulence, and is slightly different for velocity and temperature. 

This, type of scaling would not appear to correlate the measurements of either 

Morkovin or Hojstman and Owen. 

Information related to the frequency analysis of turbulent signals are 
also available from autocorrelation measurements. Autocorrelation measurements 
have been reported by Owen and Horstman, ref 47, for hypersonic boundary layers. 
Demetriades, ref. 48, has made similar measurements for supersonic wake flow. 

The autocorrelation and the energy spectra are related by a Fourier transforma- 
tion. It was also demonstrated (with the exception of verv close to the surface) 
by Owen and Horstman that the autocorrelation was approximately equal to the 
space-correlation for small times (.6 milliseconds) or convection distances 
when "Taylor's hypotheses ' 1 of space-time transformation was employed. 
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The scales of turbulence were originally defined by Taylor based on the 
space correlation coefficient, R , 

oc 

MACROSCALE L = j R dx (16) 

x 1 x 

0 

1-R 

MACROSCALE l/\ 2 = 2 LIM ( -) (17) 

X ^ ? 
x-*0 

By use of the Fourier transformation the scales are related to the wave 
number spectral function, F(k) as 

L = — F (k=o) (18) 

x 2 u 7 

and 

00 

• Jk , -L- j k 2 F ( k ) dk (19) 

X 2 u 2 ” o 

X 

Owen and Horstman, ref. 40, evaluated the turbulent macroscale from their 
autocorrelation measurements, using equation (16) and Taylor*s hypothesis. • 
Laderman and Demetriades, ref. 17 (Philco-Ford report), employed a form of 
equation (18) to evaluate the macroscale in their boundary layer. Figure 18 
is a plot of L x / 6 f° r these two hypersonic boundary layers. In figure 18a) 
the data of Owen and Horstman are for mass flux, while that of Laderman and 
Demetriades are for the longitudinal velocity. The hypersonic boundary layer 
scales are much larger, compared to the boundary layer thickness, than those 
reported for incompressible flow. Values of (L^/6) - 0.25 over the outer 
portion of the boundary layer were reported by Sandborn and Slogar, ref. 49, for 



34 


an incompressible boundary layer of 5-5 cm. The large boundary layers 
of Tieleman give values of (L^/5) - 0.1. Thus, the boundary layer thickness 
is not an adequate scaling length for the turbulent scales. 

The temperature scales measured by Owen and Horstman are smaller than 
the mass flux scales. For Laderman and Demetriades* measurements the scales 
are approximately the same. 

Evaluation of the microscale from equation (19) does not appear to have 

been reported for supersonic flows. Owen and Horstman, ref. 40, report values 

obtained by -curve fitting a parabola to the autocorrelation curves. The 

values of X^ obtained by the evaluation were nearly equal in magnitude to 

the macroscale values. From subsonic measurements the ratio X /L is 

x x 

found to be approximately 0.1 to 0.2. It is suspected that the frequency 
resolution was not adequate to define the parabola at the t = 0 intercept of 
the autocorrelation curve for Owen and Horstman's data. * 

The recent study of Owen and Horstman, ref. 40, includes data on the 
amplitude distribution of the mass flux and temperature fluctuations. Detailed 
measurements of the probability distributions of these fluctuations throughout 
the hypersonic boundary layer are presented. The mass flux fluctuation amplitude 
was found to be skewed negative near the surface, which was associated with 
large intermittent positive "spikes' 1 in the mass flux. At the outer edge 
of the boundary layer the reverse effect was found. The total temperature 
was also found to be skewed negative near the surface. Direct measurements 

_ 3/2 

of the "skewness factor" (e 1 / (e 7 ) ) are reported by Owen and Horstman. 

The intermittent aspects of hypersonic boundary layers have been measured 
by both Owen and Horstman, ref. 40, and Laderman and Demetriades, ref. 17. 
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Figure 19 shows the measured values obtained for the intermi ttency factor, y 
(where 1 is fully turbulent and 0 indicates no turbulence) , across the boundary 
layer. The intermittency factors given by Oven and Horstman was for the mass 
flow fluctuations, while those of Laderman and Demetriades also include 
temperature fluctuations. Temperature intermittency factors reported by Owen 
and Horstman were slightly higher than their mass flux values. Note also 
that the values obtained by Owen and Horstman were determined from measurements 

-rr -T 2 

of the "flatness factor," e 4 / (e 2 ) , while those of Laderman and Demetriades 
were found directly with a special electronic circuit. 

In Figure 19 it is seen that the intermittent outer region of the hypersonic 
boundary layers extend over less than 20% of the boundary-layer thickness, while 
intermittency in the subsonic boundary layer occurs over 50 to 60% of the laver. 
Both compressibility and the presence of thicker sublayers in the hypersonic 
flow may act to reduce the intermittent region of the outer flow. Reynolds . 
number is not expected to have a direct effect, as the large scale boundary 
layers of Zoric, ref. 26, and Tieleman, ref. 36, (intermittency measurements 
were made, but not reported) have intermittency factors much the same 
as measured by Klebanoff, ref. 32. Further measurements of the intermittency 
in both supersonic and high subsonic flows would be desired in order to 
determine the effect of compressibility. The early measurements of "flatness 
factor ratio" reported by Sandborn and Wisniewski, ref. 13, for a M = 3 
boundary layer were not adequate to define the intermittency factor directly 
(due to inadequate electronic multiplier response); however, it appears that 
the extent of the intermittent outer edge was quite similar to the hypersonic 
measurements shown on figure 19. 
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The supersonic and hypersonic boundary layers also contain "intermittent" 
regions near the surface. The surface intermittent region is not readily observed 
in the subsonic boundary layers. The existence of larger sublayers in supersonic 
flows, no doubt, is resoonsible for the observed intermittencv . The intermit- 
tency near the surface is also directly related to the skewness of the fluctua- 
tions found from the probability measurements. 

Space-time correlation measurements were also reported by Owen and Horstman, 
ref. 47, for a M = 7.2 boundary layer. Similar measurements in a Mach 3 turbu- 
lent wake were reported by Demetriades, ref. 48. "Convective velocities," cal- 
culated from the measured delay time, corresponding to the maximum correlation 
coefficient, were obtained from these space- time measurements. The convective 
velocities were found to be a function of the location in the boundary layer 
and of the frequency. The hypersonic convective velocity measurements appear to 
be similar to measurements made in the large-scale subsonic boundary layer, 

Cliff and Sandborn, ref. 50.* In the outer region of the boundary layer the 
convective velocities are normally less (zero pressure gradient flow) than the 
local mean velocity. For the inner region the convective velocities are found 
to be greater than the local mean velocity. The space-time correlations study 
of Owen and Horstman also give a measure of the average trajectorv of the 
turbulent fluctuations in the boundary layer. 


*[The overall convective velocities reported by Owen and Horstman are somewhat 
low due to limited frequency response of the original hot-wire instrumentation. 
Recent measurements with improved hot-wire frequency response give higher values 
for the overall convective velocity in the outer region of the flow; which are 
in better quantitative agreement with the subsonic results.] 
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d) Free Shear Flo w.- Wakes and jets are two types of free shear flows of 
major interest. The mixing of two independent flows also may be included in 
the free shear flows. While a number of studies are available for these types 
of compressible flows, they are not as well documented as the boundarv layer. 

Demetriades, ref. 48, 51 and 52, has made an extensive survey of an 
axisymmet ric wake at M =3. Values of the axial velocity, density, temperature 
and velocity-temperature fluctuations were measured. Also auto- and space-time 
correlations were evaluated. Using the proper coordinate transformation the 
compressible results were found to be "identical with their incompressible 
counterparts," Measurements, or the indirect evaluation, of the turbulent 
shear stress for a compressible wake does not appear to be available. Evaluation 
of the intermittency characteristics of supersonic and hypersonic wakes were 
reported by Demetriades, ref. 52, and by Levensteins and Krummins, ref. 53. 
Estimates of turbulent scales were also reported in these papers. ..The wake 
intermittency for compressible flow was found to be nearly the same as the 
incompressible measurements. 

A number of turbulence measurements have been reported for real gas type 
turbulent wakes; refs. 54 through 57. These data were for wakes behind 
projectiles rather than for the statistically steady flows discussed above. 
Autocorrelations, spectrum and time scales are obtained from these measurements. 
New techniques of measuring instantaneous wake profiles and turbulence are 
presented in references 58 and 59. 

Very little information on turbulence measurements in compressible jets 
was found in the present survey. Estimates of turbulent effects, such as 
eddy-viscosity or mixing lengths, are available mainly from incompressible 
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measurements. The recent paper of Peters, Chriss and Paulk, ref. 60, contain 
evaluations of the shear stress based on mean flow measurements for subsonic 
jets. Measurements of the intermittency factor for a jet with an initial 
velocity of 693 ft/sec were reported by Maestrello and McDaid, ref. 61. 

These subsonic compressible measurements agree with the incompressible inter- 
mittency data. 

A recent conference on "Free Turbulent Shear Flow," ref. 62, has examined 
in great detail existing methods of computing free shear flows. There was some 
evidence from the mean flow calculations that the compressible effects on 
turbulence may be more pronounced for the free shear flows. 

CONCLUSIONS 

A direct comparison of measured terms in the turbulent stress tensor for 
zero pressure gradient, super- and hypersonic, turbulent boundary layers with 
large Reynolds number, incompressible, boundary layers was found to show 
remarkable agreement. The turbulent stress terms were scaled with the wall 
shear stress and the boundary-layer thickness. While some variations in the 
data were found near the surface, the outer region of the zero pressure gradient, 
turbulent boundary layer appears to be surprisingly similar for all Mach 
numbers. 

Mean flow evaluations of the total shear stress distributions across the 
whole, zero pressure gradient, turbulent ,. boundary layer was also found to 
correlate nearly independent of Mach number. The total shear stress distribution 
for the compressible boundary layers was in good agreement with the measured 
Reynolds shear stress (puv) term for incompressible, zero pressure gradient. 
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layers. Recent laser anemometer measurements, by Johnson and Rose at a Mach 
number of 2.9, of puv are also in reasonable agreement with the expected 
total shear stress over the outer 60% of the boundary layer. Near the surface 
the laser measurements are systematically lower than the exnected values. 

The laser measurements are only slightly higher than the individual stress 
components (pr? and pv^") measured in incompressible flow. Thus, the direct 
laser measurements appear to be reasonable over the complete boundary layer. 

The disagreement between the measurement of puv and the total shear stress 
might be taken with some reservation to suggest that the Reynolds stress term 
may not be the only important turbulent shear stress term in supersonic 
boundary layers. 

The correlation of the measurements, together with the incompressible data 
can be proposed as a check on the super- and hypersonic measuring techniques. 
Eigure 20 is the proposed summary plot of the expected variation of the Reynolds 
stress tensor terms. Adiabatic, zero pressure gradient, turbulent, supersonic, 
boundary layers are expected to agree with these curves over the outer 50 to 
60% of the layer. Near the surface the viscous sublayer will cause deviations. 
The sublayer thickness will increase in size by greater than a factor of 10 as 
the Mach number goes from 0 to 5. 

Hot-wire anemometer measurements of the longitudinal component of velocity 
fluctuations, neglecting pressure fluctuations, for at least supersonic Mach 
numbers appear to be reasonable. However, measurements of the vertical com- 
ponent of the turbulent velocity and the Reynolds shear stress with yawed 
hot wires appear questionable. Direct measurements of the longitudinal mass 
flow fluctuations with the hot-wire anemometer show considerable variation from 
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one set of data to another. The hot wire should be most sensitive to the mass 
flow, so it appears that use of the local mean values of the mass flow are 
not adequate to normalize the results. Considerable variation is also observed 
for the temperature fluctuation measurements. 

The spectrum of mass flow fluctuation for super- and hypersonic flows 
were also found to agree with incompressible results when compared on a wave 
number plot. The resultant spectra should serve as a direct approximation of 
the frequency content of a given boundary-layer measurement. 

Measurements of the intermittency of the outer region of super- and 
hypersonic boundary layers shows a marked difference from the incompressible 
layer. The intermittency extends over a much smaller percentage of the layer 
20%) for the high-speed flows, compared to the incompressible flow (-^50%) 
Evaluation of mass-flow fluctuation intensity in the free stream of super- 
and hypersonic wind tunnels indicate the intensity increases approximately as 
Mach number squared. The level is found to decrease with Reynolds number, and 
also would be expected to decrease with increasing tunnel size. The size 
effect is not evident in the hypersonic tunnels. For Mach numbers greater 
than 2.5 the mass flow fluctuations should not be influenced by the inlet 


conditions . 
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Figure 4.- Evaluation or the zero pressure gradient, turbulent 
boundary Layer, shear stress distribution. 

M = 2.5 Meier and Rotta, ref 25 

M = 2.8 Meier and Rotta, ref 25 

M = 3.0 Meier and Rotta, ref 25 

M * 3.5 Meier and Rotta, ref 25 

M = 4.5 Meier and Rotta, ref 25 

M = 3.53 Sturek, ref 38 

M = 7.2, Rg = 4.9xi0 3 , Horstman and Owen, ref 27 

M = 7.2, R d - 7.3xl0 3 , Horstman and Owen, ref 27 

M = 7.2, Rg = 9.7xl0 3 , Horstman and Owen, new data 

The following were evaluated by Bushnell and Morris, ref 29 

M = 4.67, R e = 5.7 xIO 4 , Lee, Yanta and Leonas, ref 66 

M = 4.70, Rg = 5.0xl0 3 , Lee, Yanta and Leonas, ref 66 

M = 5.98, Rg=l.28xl0 4 , Samuels, Peterson and Adcock, ref 30 

M = 6.02, R 0 s I.40xl0 4 , Adcock, Peterson and McRee, ref 67 

M = 6.50, R e =2.53xl0 3 , Damberg, ref 28 

High Mach number data evaluated by Bushnell and Morris, ref 29 

M = 10.65 Perry, ref 68 

M = 1 2 Fiore, ref 69 

M = 20.8 Fisher, Maddalon, Weinstein and Wagner, ref 70 
M = 1.8 Gibbings and Mikulla, ref 21 
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b) Interaction region, Rose, ref 37 
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Measurement of turbulent shear stress term through a shock-wave interaction. 


Figure 7 . - 






40 X I0~ 2 


.40 X 10- 



Re 

1.52 X I0 5 
1 .89 X I0 5 
3.01 X I0 5 


V 2.8 \ 


N \ V'. 

\ \ v. 


Um E 


P u 00 


• 20 hli 


X cm Re 
.115 4900 


176 7300 


\ , 237 9700 
\ V 


% 


(a) Insulated flat plate 

Data of Meier and Rotta, Ref. 25 


(b) Cylinder with heat transfer 

Data of Horstman and Owen, Ref. 27 


Figure 9.t Turbulent heat flux evaluations for zero pressure gradient 
turbulent boundary layers. 
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a) Adiabatic 

Figure 10.- Longitudinal turbulent velocity measurements in zero 
pressure gradient boundary layers. 
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Figure 12.- Vertical turbulent velocity measurements in zero pressure 
gradient boundary layers. 
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Figure 15.-* Measurements of temperature 
gradient boundary layers. 









b) Total temperature 

Figure 18.- Macroscales of turbulence measured in hypersonic zero 
pressure gradient boundary layers. 
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Figure 19.- Measurements of the intermittency factor for zero 
pressure gradient boundary layers. 





